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We present results and prospects for searches for rare B and D meson decays with final state dimuons, in- 
cluding — -> ji + , B° — > /i + fx~ , and D° — > fi + fj,~ . Upper limits on the branching fractions are compared to 
previous CDF measurements, recent results from the B factories and theoretical expectations. We also report on 
new measurements of production and decay properties of the X(3872) particle, discovered in 2003 by the Belle 
Collaboration. New results on the measurement of the relative branching fraction for the Cabibbo suppressed 
decay B + — > J/ipn + B(B + — > J/ipn + )/B(B + — > J/ipK + ) are presented too. The presented results are based 
on the analyses of 70 to 220 pb^ 1 of data collected by the CDF II detector in pp collisions at y/s = 1.96 GeV at 
Fermilab Tevatron. 



1. INTRODUCTION 

The Tevatron collider at Fermilab continues its 
operation of Run II since 2001. The pp collisions 
are produced at y/s — 1.96 GeV with yet increas- 
ing instantaneous luminosity. The amount of in- 
tegrated luminosity delivered to the CDF and D0 
experiments has surpassed that of Run I by nearly 
a factor of 5, exceeding 0.5 fbV 1 by July 2004. 

The pp collisions are an abundant source of 
Beauty and Charm hadrons. The upgraded Col- 
lider Detector at Fermilab (CDF II) ^ can ef- 
fectively select many of the decay modes of the 
beauty and charm hadrons. This allows for an 
analysis or search of statistics limited, rare pro- 
cesses. 

2. DETECTOR 

The components of the CDF II detector per- 
tinent to the analyses presented here are de- 
scribed briefly below. Detailed description can be 
found elsewhere pQ. The cylindrical drift cham- 
ber (COT) and the silicon microstrip detector 
(SVX II) are immersed in a 1.4 T solenoidal mag- 
netic field and allow for tracking charged par- 
ticles in the range of 1 77 1 < 1 with good posi- 
tion precision (impact parameter, do, resolution 
Od ~ 15 /im) and good momentum precision 
(transverse momentum, pt, resolution Spt/p t ~ 
0.1%/(GeV/c)). The muon subdetector includes 



the CMU (detectable muon p T > 1.5 GeV/c), 
CMP (pr > 2.5 GeV/c) covering the pseudo- 
rapidity range of |7y| < 0.6, and CMX (px > 
2 GeV/c) covering the range of 0.55 < \t]\ < 1.0. 

2.1. Triggers 

CDF II employs a 3 level trigger system. Two 
sets of triggers are used by the analyses presented 
here: the dimuon and two-track triggers. 

The dimuon trigger selects two muons (track 
associated with the hits in a muon subdetector) 
with p T > 1.5 GeV/c. 

The two-track trigger selects tracks with px > 
2 GeV/c and uses the Level-2 Silicon Vertex 
Tracker (SVT) to select on the impact parame- 
ter of a track with respect to the beamline. The 
impact parameter resolution is below 50 ^m. The 
tracks are selected with 120 fim < |rf | < 1-0 mm. 

3. D° -> n+fi- SEARCH 

The flavor-changing neutral current (FCNC) 
decay D° — > is highly suppressed in the 

Standard Model (SM), B SM {D° -> p+ir) » 
10 -13 |2] ■ This prediction is many orders of mag- 
nitude beyond the reach of the present experi- 
ments. Observation of this decay at a rate signif- 
icantly exceeding the SM expectation would in- 
dicate the presence of non-SM particles or cou- 
plings. Thus a large, unexplored region exists in 
which to search for new physics. 
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This search uses a 65 pb~ data sample. The 
sample of ~ 10 5 D*-tagged two-body D° decays 
selected by the two-track trigger is used to esti- 
mate backgrounds, to optimize selection require- 
ments, and to normalize the sensitivity of the 
search from the data sample itself. 

The B(D° — > upper limit for the given 

confidence level (CL) is determined using 

B(D° — ► 7T+7T~) iV(7T7r) e(7T7r) a(7T7r) ' 

where B(D° -> tt+tt-) = (1.43 ± 0.07) x 10~ 3 
is the normalization branching fraction, N(mr) 
is the number of D° — ► ir + rr~ events observed, 
Acl^A 4 ) is the upper limit on the observed D° — ► 
signal events with the expected background 
accounted for. The e and a are the efficiency 
and acceptance for each mode. Except for the 
muon identification, and the assignment of differ- 
ent particle masses, the same selection require- 
ments are applied to both modes. Kinematically, 
the D° — > 7r + 7r~ and modes are nearly 

identical, minimizing the systematic uncertainty 
and the differences in acceptance. The efficiency 
and acceptance fractions are estimated to be 
e(7T7r)/e(/i/i) = 1.1 ± 0.04 and a(7T7r)/a(/i/i) = 
0.96 ± 0.02. The width of the reconstructed mass 
peak for two-body decays of the D° is sufficient 
to separate D° — > K~w + from ir + ir~ . 

The dominant sources of the background here 
are the combinatorial background (estimated us- 
ing the events in the high mass sideband), and 
the muon misidentification background, when a 
pion is misidentified as a muon (estimated using 
the sample of D*-tagged D° — > K~ir + events). 

A "blinded" analysis was performed. The data 
in the signal mass window were hidden and the 
analysis cuts optimized without knowledge of 
their actual impact on the result. 

Using the optimized selection requirements, 
5.0 ±2.2 events remain in the high mass sideband, 
yielding 1.6 ± 0.7 expected from the flat compo- 
nent of the background. The expected number 
of misidentification events is 0.22 ± 0.02. The to- 
tal expected background is 1.8 ± 0.7 events. The 
number of events in the normalization mode is 
N(irir) — 1412 ± 54. We apply the optimized se- 
lection requirements to the signal region of the 
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Figure 1. The mass distribution of candidate 
D° — > events. 



/i/j, sample and find no events remaining, as dis- 
played in Fig. ^ Using Eq. Q we find an upper 
limit on the branching fraction of 

B(D° -> fi+fj,-) < 2.5 x 10~ 6 (3.3 x 10~ 6 ) (2) 

at the 90% (95%) confidence level. This result 
improves on the best previously published lim- 
its. This limit has recently been improved by 
HERA-B to the value of 2 x 10~ 6 at 90% con- 
fidence level 0. We expect to improve this limit 
using a substantially larger dataset. 

4. B° (d) -> SEARCH 

The rare FCNC decay B a s n + yT is one 
of the most sensitive probes to physics beyond 
the SM 0]. The decay has not been experimen- 
tally observed and the best previously published 
branching ratio limit was B(B° — > < 
2.0 x 10~ 6 at 90% CL, while the SM prediction 
is (3.5 ± 0.9) x 10~ 9 . Similarly, the best previ- 
ously published limit on the related branching 
ratio, B(B° -> n+ >-) < 1.6 x 10~ 7 , is about 
three orders of magnitude larger than its SM ex- 
pectation. The B(B® — ► A* + aO can be enhanced 
by one to three orders of magnitude in various 
supersymmetric (SUSY) extensions of the SM. 

The presented analysis uses 171 pb _1 of data 
collected through September 2003. The data 
sample is collected using the dimuon trigger, se- 
lecting the muons with \r]\ < 0.6. 

We model the signal B°,~ — ► decays us- 

ing the Pythia Monte Carlo (MC). To normalize 
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to experimentally determined cross-section, we 
require the -B°(<f) to nave Pt(B®^) > 6 GeV/c 
and rapidity 



s(d) 

lvl<i. 



To further discriminate ^v^n 



decays 



from background events we use the following four 
variables: the invariant mass the B- 

candidate proper decay length (cr ) ; the opening 
angle (A$) between the £?-hadron flight direc- 
tion (estimated as the vector p!^ 1 ) and the di- 
rection to the decay vertex; and the B-candidate 
track isolation (I = Pj^/iYlPr^ + Pr*)' wnere 
AiZ(p§* = V(A^) 2 + (A?7) 2 < 1.0). 

We use a "blind" analysis technique. The 
data in the search mass window 5.169 < M w < 
5.469 GeV/c 2 are hidden (corresponding to w ±4 
times the mass resolution and is centered on the 
B® and B® masses) . The optimization performed 
using only data in the mass sideband regions, 
M w € [4.669, 5.169] U [5.469,5.969] GeV/c 2 . 

We vary among the sets of (M^, cr , A$, /) cri- 
teria to minimize the a priori expected 90% CL 
upper limit on the branching ratio. For a given 
number of observed events, n, and an expected 
background of n Dg , the branching ratio is deter- 
mined using: 



B(B° S -> M + M -) < 



N(n, n bg ) 
2(7^9 ^Cae 



total 



(3) 



where A^(n, «bg) is the number of candidate B® — > 
decays at 90% CL, and C is the integrated 
luminosity. The a priori expected limit is given 
by the sum over all n, weighted by the Pois- 
son probability P(n\ribg)- For the B® — > 
limit we substitute a B o for cr B o. The factor of 
two in the denominator accounts for the charge- 
conjugate i?-hadron final states. The total ac- 
ceptance times efficiency, aetotai, and rib g are es- 
timated separately for each set. 

The acceptance is estimated using the MC sam- 
ple. The total efficiency is a product of recon- 
struction, trigger, and the analysis selection. The 
reconstruction and trigger efficiencies are esti- 
mated using unbiased samples of J/ip — > 
decays. The efficiency of the analysis selection 
criteria is determined using the MC sample. 

The optimization procedure yields (cr > 
200 /j,m, A$ < 0.1 rad,J > 0.65) and the mass 
window ±80 MeV/c 2 around B° s (5.369 GeV/c 2 ) 



and B° d (5.279 GeV/c 2 ) masses. The estimated 
acceptance is a w 6.6% and the total efficiency 
is etotai ~ 30%. The expected background is 
1.1 ± 0.3 events in each of the B® and B® mass 
windows. Using the optimized set of selection 
criteria, one event survives all requirements and 
has an invariant mass of = 5.295 GeV/c 2 , 
thus falling into both the B® and B® search win- 
dows as shown in Fig. [21 Using Eq. we de- 
rive 90% (95%) CL limits of B{B° S -> < 
5.8 x 10~ 7 (7.5 x 10~ 7 ) B(B° d -► (i+u,-) < 1.5 x 
10 _7 (1.9x 10 -7 ). The new limit improves the pre- 
viously best published limit by a factor of three 
and significantly reduces the allowed parameter 
space of i?-parity violating and 5*0(10) SUSY 
models. We expect to improve the limit by us- 
ing the larger dataset as well as enhancing the 
analysis techniques [S]. 
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Figure 2. The /i + // _ invariant mass distribution 
of the events in the sideband and search regions 
satisfying all requirements. 



5. X(3872) PROPERTIES 

In 2003 the Belle Collaboration reported a new 
particle, X(3872), observed in exclusive decays of 
_B-mesons produced in e + e~ collisions 0- This 
particle has a mass of 3872 MeV/c 2 and decays 
to J/ip n + TT~ . A natural interpretation of this 
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particle would be a previously unobserved char- 
monium state, but there are no such states pre- 
dicted to lie at or near the observed mass with the 
right quantum numbers to decay into <//i/>7r + 7r _ . 
Whether it is a new form of hadronic matter or 
a conventional cc-state in conflict with theoret- 
ical models, the A (3872) is an important ob- 
ject of study. Here we report the observation 
of a J/ipTr + Tr~~ resonance produced inclusively 
in pp collisions and which is consistent with the 
A (3872). We also report the measurement of the 
A(3872) and ip(2S) production fraction associ- 
ated with a large lifetime (B-hadrons) [Jj. 

The analysis uses 220 pb _1 of data collected us- 
ing the dimuon trigger selecting the muons asso- 
ciated with CMU or CMX muon chambers within 
the mass range of 2.7 to 4.0 GeV/c 2 . 
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Figure 3. The J/?/>7r + 7r~ mass distributions for 
the sample separated into those candidates with 
dipion masses less than (circles) or greater than 
(solid points) 500 MeV/c 2 . 



pr(J/tp) > 4 GeV/c, good vertex probability, 
and AR < 0.7. 

Besides the ip(2S) peak, another peak is 
observed at the J/-i/jTr + Tt~ mass of around 
3872.3 MeV/c 2 . The Gaussian plus a quadratic 
polynomial is used to model each peak. The 
binned likelihood fit results 5790±140 ip(2S) can- 
didates and 580 ± 100 A(3872) candidates. 

The A (3872) signal reported by Belle Collabo- 
ration favors large Tr + n~ masses. Our data sup- 
ports this conclusion as well, see Fig. [3] Re- 
quiring the M„+„- > 500 MeV/c 2 reduces the 
background by almost a factor of two. We use 
the high mass sample for measuring the A (3872) 
mass, we find m(A(3872)) = 3871.3 ± 0.7(stat) ± 
OA(syst) MeV/c 2 . This is in agreement with the 
measurement by the Belle Collaboration. The ob- 
served width of 4.2 ± 0.8 MeV/c 2 is found to be 
consistent with the detector mass resolution. 

For the lifetime analysis the data sample is 
the same as for the mass measurement. The 
Af w+7r - > 500 MeV/c 2 is applied to the A(3872) 
signal region as it greatly improves the purity. 
This cut is not applied for the tp(2S) analysis, 
as it has a strong effect on the number of signal 
events due to the lower mass of the ip(2S). 
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In the offline analysis, the two opposite charge 
tracks with px > 0.4 GeV/c 2 are added to the 
J/ip candidate. All four tracks are then re- 
quired to come from the same vertex. Additional 
strict requirements are made to suppress back- 
ground: J/ip mass within 60 MeV/c 2 (^ 4cr m ), 



Figure 4. The uncorrected proper-time projec- 
tion of the full A (3872) likelihood PDF over the 
restricted mass range ±2.5 standard deviations 
around the A(3872) mass. 
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To separate the prompt from long-lived com- 
ponents we perform an unbinned likelihood fit 
in which mass and lifetime are simultaneously 
included. Each candidate is characterized by 
its mass M, uncorrected proper time, cr = 
p^(j/^) £xy, and the uncertainty on cr, <r CT . 
A projection of the fit for A(3872) is shown in 
Fig.H 

We find the long-lived fractions of the sample 
to be: F(ip(2S)) = 28.3 ± l.O(stat) ± 0.7(syst)% 
and F(X(3872)) = 16.1 ± A.%stat) ± 2.0(syst)%. 
The lifetimes of the long-lived signal components 
are consistent with the latter being produced in 
the B-hadron decays. 

6. B+ -> J/tp 7T+ MEASUREMENT 

The Cabibbo-suppressed decay B + — > 
J I ip 7r + proceeds via a b — > (cc) + d transition. 
The modes governed by this transition may show 
direct CP violating effects at the few percent 
level |B] ■ The B + — > J /ipn + is expected to have 
a rate about 5% of that of the Cabibbo-allowed 
mode B+ -> J/ip K+ . The PDG 2004 average 
of the ratio of branching ratios is 4.0 ± 0.5% [0], 
and the more recent result from BABAR is 
5.37 ± 0.37(stat) ± 0.11% 8 . We present the 
measurement of the ratio of branching ratios 
B re i = B(B+ -> J/iPir+)/B(B+ -> J/ipK+ ) 
using the CDF II data. 

Distinguishing J/tp tt + from J/tp K + events 
and determining the signal significance is a com- 
plicated task since the two modes overlap kine- 
matically. The solution is to fit in J/tpir + and 
J/tpK + mass space simultaneously. We use the 
following relation to measure the desired value: 

N(J/^n+) e J/i>K+ 1 

Wrcl — ,,, 7 , , ,,, s X — 'obs * i 

N(J/tpK+) e J/i>7 r+ £ rol 

where r Q b s is the observed raw ratio of branching 
ratios in data and e rc i is the relative efficiency of 
the two decay modes. 

The analysis uses 200 pb _1 of data collected 
using the dimuon trigger. The third track (K or 
7r) is required to have pr > 2.0 GeV/c and to be 
consistent to come from the same vertex as the 
J/ip candidate. The B + candidate should have 
Pt > 6.5 GeV/c and have a displaced decay ver- 



tex, L xy > 200 /xm. The unbinned likelihood fit 
is used to extract the sample composition. The 
J/ip 7r + and J/ip K + are each modeled with a 
Gaussian, the background in J /ip K + is mod- 
eled with a first order polynomial. The fit gives 
r ohs = 0.045±0.008with N{J/i>K+) = 1986±36 
and N(J/ip n+) = 90 ± 15. The study of the sys- 
tematic uncertainty shows that it is dominated 
by the signal model. 

The relative efficiency is measured using MC 
with full detector simulation, it is estimated to 
be e ro i = 0.991 ± 0.008. 

The measured ratio of branching ratios is 
B Iei = 4.5 ± 0.8(stat) ± 0.3(syst)%. 
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